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Three different types of maneuvers were designed to separately quantify pitch rate and 
angle of attack rate contributions to the nondimensional aerodynamic pitching moment 
coefficient. These maneuvers combined pilot inputs and automatic multisine excitations, 
and were flown with the subscale T-2 and Bat-4 airplanes using the NASA AirSTAR flight 
test facility. Stability and control derivatives, in particular C, and Cm,, were accurately 
estimated from the flight test data. These maneuvers can be performed with many types of 
aircraft, and the results can be used to increase simulation prediction fidelity and facilitate 
more accurate comparisons with wind tunnel experiments or numerical investigations. 


Nomenclature 
Ay, Az body-axis translational accelerations, g Tice body axes 
CL lift coefficient a angle of attack, rad 
Cin pitching moment coefficient B angle of sideslip, rad 
Cx, Cz body-axis force coefficients A perturbation from trim 
c mean aerodynamic chord, ft Oe elevator deflection, rad 
cov(.) covariance € downwash angle, rad 
El.) expectation operator o standard deviation 
g gravitational acceleration, ft/s? 0,9, Euler angles, rad 
I. inertia components, slug-ft? Q engine speed, rad/s 
I. engine rotor inertia, slug-ft? 
Mr thrust pitching moment, ft-lbf Subscripts 
m mass, slug 0 reference value 
N number of data points cm center of mass 
D, 9,7 body-axis angular rates, rad/s 
q dynamic pressure, Ibf/ft? Superscripts 
R? coefficient of determination - inverse 
S wing reference area, ft? - transpose 
8 Laplace variable ; time derivative 
t time, s 7 estimate 
V true airspeed, ft/s 
XT body-axis thrust component, lbf 
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I. Introduction 


HE pitch rate of an aircraft is a rate of change of the aircraft orientation with respect to an inertial 

frame, expressed in body axes. The angle of attack rate is a rate of change of the aircraft orientation 
with respect to the air-relative velocity. Arbitrary motion in flight can involve one, both, or neither of these 
rates. For relatively small perturbations, reduced frequencies, and Mach numbers, the effects of angular 
velocities on the aerodynamic forces and moments acting on the aircraft can be modeled using conventional 
stability derivatives, after Bryan. The largest contributions from these rates are on the nondimensional 
pitching moment coefficient, and are represented by the parameters C;,, and Cm, ,. The former is usually 
attributed to the additional lift acting on the horizontal tail during rotation in pitch, whereas the latter is 
an approximation to unsteady aerodynamic phenomena. The following discussion applies to other similar 
pairs of parameters, such as Cy, and Cy, but the associated results are not reported in this paper because 
the dependencies identified for these parameters were less important for the overall aircraft dynamic motion 
than the dependencies for the pitching moment. 

Model parameters Cy, and Cm, are not usually estimated separately and accurately from flight test data 
because conventional aircraft maneuvers performed for system identification analysis result in pitch rate and 
angle of attack rate time series that are nearly proportional to each other. These data are said to be highly 
correlated with each other or nearly collinear. No parameter estimation method can correctly attribute the 
dependencies of the pitching moment coefficient on these two variables when they are highly correlated. 
Solutions to this problem include combining dependencies to estimate a single combined derivative, using 
other analyses to provide prior information on some of the parameters, and redesigning the maneuver so 
that the pitch rate and angle of attack rate are less correlated. The approach taken for this paper was to 
redesign the maneuver to decorrelate the data. As discussed in the following section, Refs. 2 and 3 also used 
this approach. 

Accurate knowledge of Ci, and Cim,, is needed for prediction in flight regimes such as at high angles of 
attack, rapid maneuvering, take off, landing approaches, and loss of control scenarios. Not estimating both 
parameters in these cases can bias other model parameter estimates, which degrades the usefulness of the 
identified model. Accurate knowledge of these parameters helps to improve the fidelity of the model, and 
facilitates better comparisons with wind tunnel experiments and numerical analyses. 

This paper discusses accurate estimation of pitch rate and angle of attack rate derivatives from flight 
test data. The important characteristics of maneuvers that successfully decouple these dependencies are 
discussed, as well as some useful parameter estimation techniques. Flight test data from three different types 
of maneuvers are presented. The maneuvers are relatively simple to perform on most aircraft. Section II 
discusses the modeling and estimation problem in detail. Section III presents the experiment design and 
analysis method. The airplanes used to fly the maneuvers are described in Section IV. Flight test data and 
estimation results are presented and discussed in Section V. 


II. Problem Description 


The aerodynamic forces and moments acting on the aircraft can be modeled as truncated Taylor series 
expansions using the aircraft state and control variables for small perturbations from a reference flight 
condition. A model of the pitching moment coefficient is 

AV qe ac 
Cm = Cm + Cmy + Cm,Aat Cm, VV + Oma aV 


t Cin, Ade (1) 


where the subscript 0 denotes the reference value and A is the associated perturbation. For simplicity, 
Cm subsumes additional terms that vary with airspeed. Inherent in Eq. (1) is the assumption that the 
flow is quasi-steady, and that aerodynamic forces and moments develop instantaneously from the angles and 
velocities of the motion. This type of model is in contrast to unsteady models, e.g. Ref. 4, which depend 
upon the current and past history of the motion. 

The pitch damping derivative C,,, quantifies the pitching moment generated by a pure pitch rate while 
all other variables are constant. For conventional aircraft configurations, it is mostly attributed to the 
additional lift on the horizontal tail due to rotation about the center of mass. Etkin® derives the analytical 
expression for this effect as 


l 
C= —2a,Vin = (2) 
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where a; is the lift-curve slope of the horizontal tail, Vy is the horizontal tail volume, and I; is the distance 
between the center of mass and the aerodynamic center of the horizontal tail. 

The angle of attack rate derivative Cm, quantifies the pitching moment resulting from a sudden change 
in angle of attack, or equivalently, a pure vertical acceleration, while all other variables are constant. An 
analytical expression for this parameter similar to Eq. (2) is less straightforward because it represents an 
approximation to unsteady aerodynamic phenomena, which are more complex in nature. As discussed by 
Etkin,°® the pitching moment increment arising from an angle of attack perturbation can be written as the 
aerodynamic transfer function 


AC, (s) = (Cma +: Oms8 + OmgS? +. Fs Aa(s) (3) 


The steady-state value C,,,, Aa is reached after transients from higher-order unsteady aerodynamic terms 
have decayed. Truncating this infinite power series after the first two terms results in a quasi-steady model 
equivalent to Eq. (1). 

Another way to view the angle of attack rate derivative is the “lag of the downwash” concept by Glauert 
and Cowley.® Here it is assumed that the flow surrounding the horizontal tail is affected by the vorticity of 
the main wing following a transport delay 7 = 1,/V. The angle of attack increment at the tail is 


Aa, (t) = Aa(t) — Aa(t — 7) — (4) 


assuming a linear downwash model. Truncating the power series expansion of Aa(t — rT) in terms of 7 after 
the first two terms simplifies Eq. (4) to 


Aay(t) = (1 = 5) Aa(t) + te a(t (5) 


The expression in Eq. (4) is unsteady because it depends on previous values of angle of attack, but the 
approximation in Eq. (5) is quasi-steady because it depends on instantaneous values. Lag of the down- 
wash is an unsteady aerodynamic effect that is present in normal flight conditions, which is in contrast to 
frequency-dependent unsteady aerodynamics that are encountered at higher angles of attack.” This unsteady 
aerodynamic phenomenon is exhibited in normal flight conditions. Using Eq. (5), the analytical expression 
for the angle of attack rate derivative is 


Ome = —20:Ve = — (6) 


While Eqs. (2) and (6) differ only by the downwash derivative, the physical origins from which pitch damping 
and the angle of attack rate derivative arise are different. The angle of attack rate derivative can also absorb 
or mask effects due to aeroelasticity, compressibility effects, and changes in center of mass. 

The accurate estimation of model parameters in Eq. (1) from flight test data is a common goal for 
system identification analyses. Conventional maneuvers performed to obtain flight test data involve small 
perturbations about a reference condition, e.g. using frequency sweep, multisine, or step-type inputs. For 
these types of maneuvers, pitch rate and angle of attack rate time series are almost identical to each other. 
As discussed in Ref. 2, this correlation can be understood from examining the differential equation for the 
angle of attack 


Xr sina 
mV cos 8 


qs _ g (cosa cos ¢ cos @ + sin asin @) 
mvV cos 6 es V cos 6 


tan B (pcosa +rsina) — (7) 
The terms on the right-hand side of Eq. (7) are ordered in descending magnitude. The first term is the 
pitch rate, arising from kinematic terms in the equations of motion. The second and third terms are the 
contributions from aerodynamic lift and gravity forces. The last two terms arise from lateral kinematics and 
thrust variations, both of which are typically small. For small perturbation maneuvers, the kinematic pitch 
rate term dominates the angle of attack rate equation. It is therefore expected that g and @ are almost equal 
and therefore highly correlated in these types of maneuvers. 

To illustrate this correlation, flight test data for a typical small perturbation maneuver using the T-2 
aircraft (described in Section IV) is shown in Fig. 1. Multisine inputs (described in Section III) were applied 
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to the elevator, aileron, and rudder surfaces while the aircraft was in a straight and level flight condition to 
elicit small perturbation responses about all axes of the aircraft. The input bandwidth was 0.2 to 2.2 Hz, 
which spanned the expected short period and Dutch roll frequencies. This maneuver produced good ranges 
and signal-to-noise ratios in the response variables for linear modeling. Angle of attack rate was computed 
by correcting wing-tip vane measurements to the center of mass and smoothly differentiating® the data. A 
visual inspection of gq and a in Fig. 1 shows they look similar. This notion is confirmed by the plot matrix in 
Fig. 2. The diagonal of this matrix labels the explanatory variables corresponding to the model parameters. 
The upper triangle consists of pairwise cross-plots of the data. Units are not indicated because the shape of 
the data is examined rather than the magnitude. Most combinations of variables cover the space well and do 
not reflect a strong linear dependence. However, the pitch rate and angle of attack rate move in proportion 
to each other and their cross plot resembles a thin straight line, which visually indicates their correlation. 
Pairwise regressor correlations for two series, x; and x, are defined as 


Tjk = = (8) 
>= (aj () — 23)" ve (we (i) — Bx)? 


where £ denotes the mean of regressor x. This metric ranges between +1 and is a measure of linear depen- 
dence between the pair of measurements. Correlations with absolute value greater than 0.9 are considered 
highly correlated. Figure 2 shows that only the explanatory variables corresponding to Cm, and Cm, ex- 
hibit high correlation. Diagnostic tests*:? indicated that multicollinearity was not an issue, so that pairwise 
correlations are sufficient in diagnosing the collinearity. 

Figure 3 shows time histories of the angle of attack rate and the first three terms on the right-hand side 
of Eq. (7). The lift coefficient was computed as 


Cee ae as. (9) 
Maz 

Cz = 10 

z= Fs (10) 

Cr = —Cz cosa + Cx sina (11) 


which used measured data and is independent of an assumed aerodynamic model structure. The kinematic 
pitch rate contribution is dominant and describes 83% of the total angle of attack rate. The contributions 
from the lift force make a small difference, and the contribution from gravity is almost constant. 

Correlated data is problematic for estimating model parameters because variation in the dependent 
variable cannot be accurately attributed to the explanatory variables. The resulting parameter estimates 
will be biased and have increased uncertainties. One solution to this problem is to supply prior information to 
the estimator from a separate analysis, such as wind tunnel tests, using a Bayesian formulation.® Care must 
be taken, though, because any error in the supplied information diminishes the accuracy of the remaining 
estimated parameters. A second solution, which is routinely used, is to adjust the model structure so that 
the explanatory variables are not correlated. For the pitching moment coefficient, the pitch rate and angle 
of attack rate dependencies are combined into the one parameter, Cy, + Cm,- The total variation is then 
attributed to pitch rate, rather than angle of attack rate, for ease of implementation in state-space models. 
Attributing this variation instead to the angle of attack rate would result in an implicit equation in a. 
Lumping the dependencies in this way will bias parameter estimates and degrade predictive capability for 
flight regimes where pitch rate and angle of attack rate are no longer correlated. Furthermore, models 
based on numerical or analytical predictions cannot be reliably updated from flight test results because the 
parameters are confounded. A third solution, which is the approach used in this paper, is to redesign the 
maneuvers such that the data are less correlated. 

Continuing the T-2 example, Fig. 4 shows the parameter estimation results from five repeated small- 
perturbation maneuvers, with pitch rate and angle of attack rate dependencies estimated separately and as 
a lumped parameter. Results are shown using the equation-error method described in Section III, and error 
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Figure 1. Measurement time histories, T-2 straight and level maneuver 


EE, 
0.02 Aa Pay 
0.09 0.21 
0.02 0.00 0.92 
0.00 ~0.20 ~0.83 0.70 Abe 


Figure 2. Regressor cross plots and pairwise correlations, T-2 straight and level maneuver 


5 of 20 


American Institute of Aeronautics and Astronautics 


— Total a 
- Pitch rate 
30 — Lift 
» Gravity 


20 


deg/s 0 


Time, s 


Figure 3. Angle of attack rate contribution time histories, T-2 straight and level maneuver 


bars reflect the 2o uncertainty bound. The variation in airspeed was small, so Cy, was not estimated. In 
the first case, Cy, was not included in the model structure and C,,, subsumed both pitch rate and angle 
of attack rate dependencies. The parameter estimate scatter and uncertainty bounds were both low. In the 
second case, C;,, was included in the model structure. The estimates of C,,, and Cm,. were only slightly 
affected by whether or not C,,, was included in the model. However, the estimates of Cm, incurred a 
significant bias, and estimates of both C,, and Cm, had large uncertainties and large amounts of scatter 
between maneuvers. 
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Figure 4. Parameter estimates, T-2 straight and level maneuvers 


The nonlinearities present in Eq. (7) offer creative opportunities for designing maneuvers that reduce 
the amount of correlation between angle of attack rate and pitch rate. Maine and Iliff? investigated two 
maneuvers that used large bank angles to decorrelate pitch rate and angle of attack rate through the gravity 
term in Eq. (7). The first maneuver, called a smooth roll, was a steady 360 deg roll performed between 
30-60 deg/s while elevator pulses were applied. The second maneuver, called a two-point hesitation roll, 
involved applying an elevator pulse with wings level, then rolling the aircraft 180 deg and repeating the 
elevator pulse while inverted. In both these maneuvers, cos @ varied between 1 and —1, which decorrelated 
pitch rate and angle of attack rate. Parameter estimation results were repeatable and in rough agreement 
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with simple analytical predictions. It was noted that this maneuver was performed using a small and agile 
jet trainer, and that it may not be acceptable for other aircraft. 

Jategaonkar and Gopalratnam? presented an alternative maneuver that also exploited the gravity term in 
Eq. (7) but avoided rolling the aircraft into an inverted position. The maneuver involved applying multistep 
elevator inputs while in straight and level flight, and then also while banked to at least 40 deg. The throttle 
was also moved to vary airspeed. These changes were sufficient to reduce the correlation between q and a, 
and estimated values of the model parameters were similar to wind tunnel predictions. 

Any maneuver that sufficiently decorrelates gq and @ can be used for identification. Successful maneuvers 
in the past have relied upon achieving a large bank angle and using gravity to achieve the decorrelation, as 
can be understood from Eq (7). In this paper, other strategies and maneuvers are used to decorrelate the 
pitch rate and angle of attack rate. The maneuvers are applicable for testing with most aircraft. Practical 
aspects of the equation-error method are discussed, although any other parameter estimation method could 
be used. The focus is on achieving good flight test data for modeling. 


III. Methods 


This section briefly summarizes a few pertinent test and analysis techniques. Software for the input 
design, post-flight data compatibility analysis, signal processing, data analysis, and parameter estimation 
are available in a MATLAB® toolbox called System IDentification Programs for AirCraft, or SIDPAC,!° 
which is associated with Ref. 8. 


III.A. Orthogonal Phase-Optimized Multisine Excitation 


Some maneuvers discussed used orthogonal phase-optimized multisines, which are perturbation inputs 
developed at NASA Langley Research Center to achieve efficient excitation of multiple inputs for system 
identification analysis. References 8, 11, and 12 provide details on the input design, and Ref. 13 presents 
numerous flight test examples. 

Each perturbation input Au;, to be applied to input u,;, is a summation of sinusoids of the form 


ke Ky 


The selection of the maneuver duration T determines the fundamental frequency of the excitation, 1/T. 
Multiples k, chosen over a desired frequency bandwidth, then define the excitation frequencies 2rk/T. 
Because harmonic sinusoids with different frequencies are orthogonal in both the time and frequency domains, 
the harmonics k can be assigned to different inputs and applied simultaneously without correlating the inputs. 
The relative power at each frequency can be chosen arbitrarily using a,. Once these parameters have been 
designed, the phase angles ¢, are optimized to produce inputs with minimal relative peak factors, to keep 
the aircraft relatively close to the desired flight condition. The input is scaled using gain A; according to 
desired response levels and signal-to-noise ratios. 


IlI.B. Equation Error Parameter Estimation 


Given adequate data, many different methods of parameter estimation can accurately identify the model 
parameters in Eq. (1). For simplicity, equation error and time-domain data were used in this work. The 
method is briefly introduced here, and more details can be found in Ref. 8. Equation error is often reduced 
to an ordinary least squares problem, i.e. 

z= X0+Vv (13) 


where z is the dependent variable, X is a matrix of regressors or explanatory variables, 0 is a vector of 
unknown constants to be estimated, and v is a vector of equation errors. The cost function 


J(0) = ; (z — X0)" (2 — X86) (14) 


is minimized by the optimal solution for the model parameters 


6 = (X7X) x7: (15) 


7 of 20 


American Institute of Aeronautics and Astronautics 


Equation (1) can be rearranged into the canonical form of Eq. (13) for estimating the model parameters. 
The dependent variable is the nondimensional pitching moment coefficient, which can be computed from 
measured data as 


1 


Tyyq + Ulex — Izz) pr + Inz (p? — 1?) — IeQer — Mr] (16) 
The last two terms account for rotational momenta and thrust moments from rotating engines. The pitch 
acceleration g was numerically computed using a smoothed differentiation technique.® Equation (16) accounts 
for the full nonlinear dynamics of the aircraft motion, as well as time-varying dynamic pressure and mass 
properties. The regressor matrix and model parameter vector are 


1 SYD vail) Se HS Aee(1) . 
AV (2) (2)e a(2)E Cm 
oe 1 =j> Aal2) Hay ayy Ade(2) 6=|c " (17) 
AVN) a(N)e A(N)e | 
1 =j- Aa(N) syny svny Ade(N) Cms. 


For the T-2 flight test data in straight and level flight, the bottom plot in Fig. 1 shows the calculated pitching 
moment and the model output computed using equation error. The match was very close and had coefficient 
of determination R? = 0.95, indicating a good fit to the data. 

Etkin® discusses normalizing the rate variables by ¢/2Vo in the small disturbance case, and by ¢/2V 
in the general case. Using the reference airspeed Vo in the normalization can reduce correlation of the 
other explanatory variables with airspeed for small perturbation maneuvers. However, as shown in Fig. 2, 
correlations with airspeed are among the lowest in small perturbation maneuvers. Because some of the 
maneuvers discussed involve appreciable changes in airspeed, ¢/2V was used for normalization. No noticeable 
differences were observed in the results for the small-perturbation maneuvers. 

Ordinary least squares theory assumes that the regressor matrix X is measured without error. For this 
reason, an effort was made to improve the data prior to estimation. A data compatibility analysis® was 
employed to remove time skew, bias, and scale factor errors in the measurements. The regressors were also 
smoothed to discard measurement noise, which would bias parameter estimates. Smoothing the data also 
helps to illuminate and diagnose multicollinearity. 

The uncertainty in the parameter estimates is* !4 


cov(6) = (XTX) | XT E[vv]X (XTX) * (18) 


where the sample autocorrelation function is used to estimate the covariance of the error from the residuals. 
This formulation is more complex than that typically used in ordinary least squares, and arises from serial 
correlations in the data due to model structure error, which is not assumed by the theory. 

Sometimes there is not enough information content in a single maneuver to estimate all the model 
parameters of interest. When using equation error, data from multiple maneuvers can be combined easily, 
and in multiple ways,® to increase the data information content. One simple procedure is to concatenate the 
data from the maneuvers. For m maneuvers, Eq. (13) is augmented to become 


Zi XxX, Vi 
Z xX Vv 

ela] oF jot] 2 (19) 
Zm Xm Vn 


and the analysis then proceeds in the usual fashion. Combining data from multiple maneuvers is one way 
to account for data information content deficiencies. This technique assumes that the model structure is the 
same, and that the error variances are similar, for all maneuvers. When the error variances differ among 
the maneuvers, a weighted average results. This could be the case when different types of maneuvers are 
combined, because the error variance includes model structure error, which is maneuver dependent. 
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IV. Aircraft and Flight Tests 


The T-2 and Bat-4 aircraft were used to fly the maneuvers. These airplanes are pictured in Fig. 5 with 
mass and geometry properties given in Table 1. Both aircraft were tested using the NASA Langley AirSTAR 
(Airborne Subscale Transport Aircraft Research) flight test facility, where a research pilot flew the aircraft 
remotely using synthetic vision drawn from telemetry data and a map of the local terrain. Multisine inputs 
were applied by the pilot pushing and holding a button on the thrust lever while flying the airplane. Data 
were telemetered at 200 Hz, and downsampled to 50 Hz for real-time modeling. For more information on 
AirSTAR, see Ref. 15. 


(a) T-2 airplane (credit: NASA Langley Research Center) (b) Bat-4 airplane (credit: NASA Wallops Flight Facility) 


Figure 5. Test aircraft 


Table 1. Test airplane geometry and nominal mass properties 


Parameter Unit T-2 Bat-4 


b ft 6.849 12.68 
é ft 0.915 1.500 
S tf? 5.902 19.02 
ces in 57.30 48.98 
ee, in 0.000 0.036 
oa in 11.43 24.01 
m slug 1.585 3.593 


slug-ft?2 1.179 9.841 
slug-ft? 4.520 8.861 
slug-ft? 5.527 16.17 
slug-ft? 0.211 —1.015 


8 


Xx 


ars 


Xx 


The T-2 aircraft is a 5.5% dynamically-scaled model of a generic, commercial, jet transport aircraft. 
It has twin jet engines mounted under the wings and retractable tricycle landing gear. The airplane was 
equipped with a micro Inertial Navigation System (INS), which provided 3-axis translational accelerometer 
measurements, angular rate measurements, estimated attitude angles, and Global Positioning System (GPS) 
velocity and position. An inertial measurement unit (IMU) gave additional accelerometer and gyroscope 
measurements with lower latency. Air data probes attached to booms mounted on each wingtip measured 
angle of attack, flank angle, static pressure, and dynamic pressure. Measurements from static pressure 
sensors and ambient temperature sensors were used to compute air density and altitude. Engine speed was 
measured using a tachometer and used as input to an engine model to compute thrust. The engine model was 
identified from ground test data with adjustments for ram drag identified from flight data. Potentiometers 
on the rotation axes of the control surfaces were used to measure control surface deflections. Mass properties 
were computed based on measured fuel flow, pre-flight weight and balance, and inertia measurements done 
on the ground for the aircraft without fuel. The T-2 aircraft has 16 separate control surfaces. For the flight 
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data analyzed in this work, only the elevator, aileron, and rudder control surfaces were deflected. 

The Bat-4 is a fixed-wing, remotely-piloted, unmanned aircraft. It has fixed tricycle landing gear, a 
twin-boom tail section, a single fixed-pitch pusher propeller, and a 6.6 in® two-stroke gasoline engine. It is 
manufactured by the MLB Company and marketed as an aerial mapping, surveillance, and payload plat- 
form. The airplane was fitted with custom hardware and sensors for research flight testing. Onboard sensors 
included rate gyros, translational accelerometers, outside air temperature sensors, GPS, and estimated Euler 
angles. Vanes on dual wing-tip mounted air-data booms provided angle of attack and flank angle mea- 
surements. Pressure taps on the booms also measured static and total pressures, from which airspeed and 
altitude were computed. Fuel flow measurements were combined with a CAD model of the aircraft and pre- 
flight weight and balance data to model the changing mass properties during flight. The ailerons were moved 
asymmetrically to apply aerodynamic roll moments, and the ruddervators were moved both symmetrically 
and asymmetrically to apply aerodynamic pitch and yaw moments, respectively. 


V. Results and Discussion 


V.A.  q-Isolation Maneuvers 


Traditional pullup or pushover maneuvers involve the pilot slowly pulling back or pushing forward on 
the stick. Similar pullup maneuvers and pushover maneuvers were flown to separate pitch rate and angle of 
attack rate, with the important difference that angle of attack was held approximately constant during the 
maneuver. This was done to create a zero angle of attack rate with a non-zero pitch rate, consistent with 
the definition of the C,,, parameter. 

These qg-isolation maneuvers were developed using the AirSTAR piloted flight simulator, and were refined 
and practiced during several simulation sessions prior to the flight tests. The most time-consuming part of 
the maneuver development was determining the initial conditions for the entry into the maneuver and for the 
start of the constant angle-of-attack phase. Once the maneuvers were developed and refined in the simulator, 
the two flight test maneuvers were completed successfully, with one repeat each, and provided excellent data 
for modeling. Having an angle-of-attack indicator at the pilot station and the use of an angle-of-attack- 
command flight control law made the piloting task relatively easy. It is believed that the piloting task would 
also be manageable with other flight control laws, including direct stick-to-surface control in un-augmented 
aircraft, as long as the pilot has a real-time display of angle-of attack at the pilot station. 

For the pullup maneuver, the airplane was initially trimmed in level flight at an indicated airspeed of 70 
knots. Power was set to 50% RPM, the nose was smoothly pushed over to a target pitch attitude of —20 
deg, and then back stick was applied to capture and hold 6.0 deg angle of attack. This angle of attack was 
held as the airplane slowly pitched up, and recovery was initiated at a pitch attitude of 25 deg. 

For the pushover maneuver, the airplane was initially trimmed in level flight at an indicated airspeed of 
90 knots. Power was set to 66% RPM, the nose was smoothly pulled up to a target pitch attitude of 25 deg 
nose up, and then the stick was relaxed to capture and hold 4.5 deg angle of attack. This angle of attack 
was held as the airplane slowly pitched down, and recovery was initiated at a pitch attitude of —20 deg. 

Figure 6 shows flight test data for q-isolation pullup and pushover maneuvers. Time t = 0 in these plots 
corresponds to when the target angle of attack was captured. These plots show that a pitch rate registered 
while the other quantities were held approximately constant or varied slowly. The contributions to the angle 
of attack rate are shown in Fig. 7. The total value is almost zero throughout the maneuver because the 
angle of attack was mostly constant. The other contributions cancel out when summed together, and the 
kinematic pitch rate was less dominant than in the small-perturbation maneuver. 

By design, pitch rate and angle of attack rate had a low level of correlation during this maneuver. This 
is confirmed by the cross plots in Fig. 8 and the low pairwise correlations between Cm, and Ci,,, of 0.15 and 
0.10. However, these plots also indicate that the airspeed, pitch rate, and elevator were correlated. This is 
because all of those quantities changed slowly and in a similar manner. These q-isolation maneuvers could 
therefore not be used to independently estimate all the model parameters in Eq. (1). 

However, because the correlations in the q-isolation maneuvers are complementary to those in the straight 
and level maneuvers, the two maneuvers were combined in the estimation analysis to achieve independent 
estimates of all the model parameters. The straight and level maneuvers provided information on the 
angle of attack, elevator, and combined pitch rate and angle of attack rate. The gq-isolation maneuvers 
provided information on the pitch rate, with angle of attack rate held approximately constant. Between the 
two maneuvers, all model parameters were estimated accurately. The fit to C,, shown in Fig. 6 includes the 
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Figure 7. Angle of attack rate decomposition, T-2 qg-isolation maneuvers 
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Figure 8. Regressor cross plots and correlation coefficients, T-2 q-isolation maneuvers 


perturbation maneuvers at straight and level flight conditions. Analyzing different maneuvers simultaneously 
is similar to the approach used in Ref. 2. The difference here is in the new maneuver that is used to decorrelate 
the explanatory variables. 

Figure 9 shows the estimated model parameters and associated 20 bounds. The first two maneuvers are 
q-isolation pullups and the last two are q-isolation pushovers. Each result was obtained using one q-isolation 
maneuver and five straight and level maneuvers. Fewer straight and level maneuvers could have been used, 
but this was simple to implement and improved the parameter uncertainties without significantly altering 
the parameter estimates. Estimates using all nine maneuvers together are given in Table 2, and are similar 
to those illustrated. The identified parameters show good consistency between maneuvers, and relatively 
small error bounds for most parameters. Some of the scatter in the C,,,, estimate is because this term 
subsumes the quantity Cy, +2C,, and C,,, is different between maneuvers. From evaluating Eqs. (2) 
and (6), values of —53.8 and —10.7 rad~! were expected for the pitch damping and angle of attack rate 
derivative. These values are similar to and in statistical agreement with the identified parameters. Wind 
tunnel testing using a similar aircraft model reported that unsteady aerodynamics for this type of aircraft 
at low angles of attack are negligible.!° Therefore C,,,, should be attributed to lag of the downwash, which 
corroborates these results. 


Table 2. Parameter estimates and uncertainties 


Parameter T-2 Bat-4 
Cmy +0.05+0.01  —0.01 + 0.00 
Cm —1.50+0.03 —0.23+0.01 
Cm —52.643.24 —7.96+0.58 
Cme —7174+2.23 +2.67+0.35 
Cms —1.79+0.04 —0.28+0.01 


e 


V.B. Multiple-Mode Excitation 


The response of the classical airplane phugoid mode is primarily characterized by a slow oscillation in 
airspeed, pitch rate, and pitch angle at approximately constant angle of attack. Exciting this modal response 
in flight decorrelates pitch rate from angle of attack rate because the angle of attack remains relatively 
constant, similar to the q-isolation maneuver. For estimating both pitch rate and angle of attack rate 
dependencies, however, this type of maneuver would normally have to be combined with another maneuver, 
also similar to the q-isolation maneuvers, to gain information for estimating all the model parameters. 

Multisine inputs were designed for the Bat-4 to move the ruddervator and aileron surfaces in different 
ways. The inputs had a 40 s duration, uniform power spectra, and contained frequencies between 0.05 and 
1.5 Hz, which excited both phugoid and short period responses. These inputs were applied after the aircraft 
had been trimmed in straight and level flight. The longer maneuver duration and lower excitation frequencies 
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Figure 9. Parameter estimates, T-2 straight and level and q-isolation maneuvers 


were possible because this was the first flight test campaign in which the AirSTAR ground control station 
was used to fly an unmanned vehicle beyond visual range of the safety pilot.!” 

Flight test data for this maneuver are shown in Fig. 10. The elevator is considered the symmetric motion 
of the two ruddervator surfaces. The lower frequencies of the multisine input excited the phugoid mode, 
which is primarily evident in the airspeed and pitch angle measurements, and the higher frequencies excited 
the short period mode, which is apparent in the angle of attack and pitch rate measurements. The inputs 
created small perturbations from the reference condition. 

Figure 11 shows the regressor cross plots and pairwise correlations. The experiment design resulted in 
enough phugoid motion so that pitch rate and angle of attack rate were decorrelated to a workable extent, 
indicated by a pairwise correlation coefficient of 0.56. All other correlations were low to moderate, and the 
explanatory variables covered good ranges for modeling. Diagnostic metrics did not indicate a concern for 
multicollinearity. 

In Fig. 12, the angle of attack rate is decomposed into its component parts. Although the kinematic pitch 
rate was still the most dominant term, changes in the lift and gravity terms were large enough to lower the 
correlation between pitch rate and angle of attack rate sufficiently for accurate parameter estimation. The 
gravitational contribution changed because the airspeed and pitch angle, in addition to the other variables, 
varied. Likewise, the lift contribution changed because the lift force and airspeed varied. Although only 
about two cycles of the phugoid mode were evident in the maneuver, this duration was long enough to 
sufficiently reduce the regressor correlations. 

Parameter estimates are shown in Fig. 13. In the first three maneuvers, each aileron and each ruddervator 
control surface was moved independently. In the fourth maneuver, the control surfaces were moved in the 
traditional way. Different maximum amplitudes on the deflections were applied in each maneuver. The 
parameter estimates are in statistical agreement among the four maneuvers. Based on Eq. (2), Cm, was 
expected to be about —11.8 rad~!. This approximation is questionable because of the unusual airplane 
design, but is in rough agreement with the results shown in Fig. 13. The estimates for C,,, were positive, 
which is not consistent with the lag of the downwash concept. The physical mechanism for this result 
is currently not understood, but interactions resulting from the pusher propeller and tail being in close 
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proximity are suspected. Regardless of the parameter estimates for this specific airplane, the maneuver 
resulted in flight test data that had low correlations and good quality for modeling. Parameter estimates 
using all four maneuvers together are given in Table 2. 
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Figure 10. Measurement time histories, Bat-4 phugoid excitation maneuver 
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Figure 11. Regressor cross plots and correlation coefficients, Bat-4 phugoid excitation maneuver 


V.C. Excitation Through Turns 


Maine and Iliff? used a 360 deg roll to separate pitch rate and angle of attack rate. Jategaonkar and 
Gopalratnam? also used this strategy, but reduced the bank angle to 40 deg and added a segment of data in 
straight and level flight. Similar information can be extracted from data obtained from flight through a turn. 
Before AirSTAR was equipped with beyond visual range capability, operation in line of sight of the safety 
pilot restricted the aircraft to low maximum altitudes and ranges. Flight paths therefore often consisted of 
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Figure 12. Angle of attack rate decomposition, Bat-4 phugoid excitation maneuver 
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Figure 13. Parameter estimates, Bat-4 phugoid excitation maneuver 
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a short straight section followed by a 180 deg turn. 

In the maneuvers that follow, multisine inputs were applied to the T-2 aircraft either continuously or 
intermittently as the aircraft flew a trajectory similar to a typical oval-shaped holding pattern. Flight data 
for one maneuver is shown in Fig. 14. The first 15 s consist of flying a straightaway with multisine excitation. 
The next 25 s of data were in a 45 deg bank angle turn, with the multisine still running. Normal operational 
limits restricted the bank angle from increasing much further. To maintain altitude, the pilot pulled back 
on the stick and increased angle of attack. 

The correlation plots in Fig. 16 show that pitch rate and angle of attack rate had a correlation of 
0.64, which was sufficiently low for accurate parameter estimation. All other correlations were also low to 
moderate, and there was no indication of multicollinearity in the data. Figure 15 shows that pitch rate and 
angle of attack rate became most different between 15 and 40 s, when the aircraft was banked to 45 deg 
and the gravitational contribution to angle of attack rate changed the most. Although this maneuver was 
sufficient for modeling, larger bank angles would have further decorrelated these quantities. 

Parameter estimates are shown in Fig. 17. The first two maneuvers are continuous multisine excitation 
through the turns, and with moderate turbulence. The second two maneuvers are intermittent excitation 
through the turns, and with light turbulence. Parameter estimates vary somewhat between maneuvers, 
particularly Cy,,,, because the angle of attack ranged between 2 to 16 deg. This large range in angle of 
attack includes known nonlinear phenomena such as stall and a nose-slice departure, and abuses the linear 
model structure assumed in Eq. (1). However, estimates are similar to those shown in Fig. 9, and are in 
statistical agreement with analytical predictions. 
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Figure 14. Measurement time histories, T-2 turn excitation maneuver 
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Figure 15. Angle of attack rate decomposition, T-2 turn excitation maneuver 
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Figure 17. Parameter estimates, T-2 turn excitation maneuver 


18 of 20 


American Institute of Aeronautics and Astronautics 


VI. Concluding Remarks 


This paper discussed maneuvers suitable for separately quantifying the effects of pitch rate and angle 
of attack rate on the nondimensional aerodynamic pitching moment coefficient. Three types of maneuvers 
were presented. The first type was a q-isolation pullup or pushover maneuver flown at constant angle 
of attack. These maneuvers had nonzero pitch rate but approximately zero angle of attack rate. When 
combined with other maneuvers flown in straight and level flight, excellent modeling results were obtained. 
The second maneuver type involved exciting both the phugoid and short period modes using an optimized 
multisine input. The low-frequency excitation resulted in an airspeed variation that contributed towards 
decoupling pitch rate and angle of attack rate. The higher frequency excitations provided the remainder of 
the data information content necessary for estimating all model parameters. The last maneuver type applied 
optimized multisine excitations in both straight and level flight and in turning flight. The 45 deg bank angle 
was sufficient in decorrelating the data for modeling. 

The maneuvers presented were relatively easy, if not routine in some instances, to accomplish. They 
can be performed with almost any aircraft instrumented for applying stored excitation inputs on the control 
surface commands, and did not require achieving large attitude angles. The maneuvers discussed are only 
a few of the possible maneuvers that can be used to decorrelate the data. Equation (7) suggests several 
possibilities for flying maneuvers that decorrelate the data. 

The problem addressed here is a data information content deficiency; it is not a problem with parameter 
estimation per se. Past work has relied on the using the output-error method for parameter estimation. This 
work used the equation-error method, which is simple, fast, and provides a range of statistical metrics to the 
analyst. It also provides a simple means for combining maneuvers in the analysis, which was useful for the 
q-isolation maneuvers. 

The new maneuvers presented, combined with the suggested parameter estimation techniques, represent 
a way to accurately separate pitching moment dependencies on pitch rate and angle of attack rate in an 
efficient manner. This information can be used to more precisely compare flight test results with other 
results, such as from simple analytical predictions, numerical aerodynamic codes, and wind tunnel tests. 
These maneuvers can also aide in building realistic flight simulation databases for unusual flight conditions 
where pitch rate and angle of attack rate affect the aircraft motion in different ways. 
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